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ABSTRACT 


This  report  describes  the  mathematical  models  and  algorithms  that 
comprise  the  ballistic  tristatic  tracking  filter  used  by  the 
Multistatic  Measurement  System  at  the  Kiernan  Re-entry 
Measurement  Site. 
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INTRODUCTION 


This  report  describes  the  mathematical  models  and  algorithms 
that  comprise  the  tristatic  tracking  filter  used  by  the 
Multistatic  Measurement  System  (MMS)  at  Kiernan  Re-entry 
Measurement  Site  (KREMS).  The  rest  of  this  section  gives  a  brief 
description  of  MMS.  Section  2  discusses  the  basic  physics  and 
mathematical  model  upon  which  the  tristatic  filter  is  based. 
Section  3  details  the  implementation  of  this  Kalman  filter. 
Section  4  compares  the  tracking  accuracy  of  the  tristatic  filter 
to  the  accuracy  of  the  TRADEX  monostatic  tracker.  Some 
mathematical  details  have  been  collected  in  the  appendixes. 

1.1.  Basic  Concept  of  the  Multistatic  Measurement  System 

The  Multistatic  Measurement  System  consists  of  the 
monostatic  radar  TRADEX  and  two  remote  sites  on  the  islands  of 
Gellinam  and  Illeginni. 

A  microwave  link  provides  real  time  communication  between 
the  three  sites  via  a  repeater  station  at  Gagan.  See  Fig.  1. 

The  equipment  at  each  remote  site  includes: 

(1)  A  radar  antenna  and  receiver. 

(2)  Microwave  communications  hardware. 

(3)  A  clock  used  to  time  tag  the  radar  returns. 
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Configuration  of  Multistatic  Measurement  System 


Control  signals  are  sent  between  TRADEX  and  the  remote  sites 
via  the  microwave  link  to  synchronize  the  remote  site  clocks  with 
the  one  at  TRADEX.  When  TRADEX  transmits  a  radar  pulse  the 
signal  reflects  from  the  target  to  the  TRADEX  receiver  and  to 
each  remote  site  receiver.  The  propagation  times  from  the  TRADEX 
transmitter  to  each  of  the  three  receivers  measure  the  TRADEX 
range  RT,  the  Gellinam  bistatic  range  sum  R-j+Rq,  and  the 
Illeginni  bistatic  range  sum  R-p+Rp .  The  remote  sites  use  the 
microwave  link  to  report  these  measurements  to  the  KREMS  Data 
Center  ( KDC)  in  real  time.  TRADEX  also  reports  its  range  R-p, 
azimuth  and  elevation  to  KDC  in  real  time. 

The  KDC  computer  provides  two  major  functions  for  MMS.  It 
controls  the  equipment  at  the  remote  sites  via  the  microwave 
link,  and  it  executes  the  tristatic  filter.  The  results  from 
this  filter  can  be  used  to  direct  the  remote  sites. 
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GENERAL  MATHEMATICAL  MODEL  OF  MMS 


This  section  presents  the  mathematical  model  used  in  the  V 

filter  for  the  target  motion  between  measurements  and  the  7 

formulae  relating  target  position  to  the  MMS  measurements. 

i*v 

2.1.  Coordinate  System 

The  coordinate  system  used  in  this  filter  is  the  KDC 

reference  system.  This  is  a  cartesian  system  with  origin  at  the 

:  k- 

TRADEX  antenna  oriented  so  that  the  x,  y,  and  z  axes  are  pointing  -  ; 

eastward,  northward,  and  upward  respectively.  This  choice  of 

coordinate  system  allows  for  easy  computation  of  ranges  from  the  - 

L. 

three  antennas  using  their  coordinates.  It  is  also  a  convenient 
coordinate  system  in  which  to  express  target  motion. 

2.2.  Equations  of  Motion  ^ 

L 

The  equations  of  motion  discussed  in  this  section  are  used 
in  the  filter  to  predict  the  motion  of  the  target  between  7 

measurements.  + 

i 

The  equations  of  motion  can  be  written  as  a  system  of  first 
order  differential  equations  using  the  state  vector  7 

x  =  (x,  y,  z,  x,  y,  z)T.  In  vector  form,  these  equations  are 

i 
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The  first  three  components  of  this  equation  merely  state 
that  the  derivative  of  target  position  is  its  velocity.  The 
bottom  three  components  state  that  the  derivative  of  velocity  is 
acceleration. 

This  acceleration  is  expressed  in  three  parts.  The 
acceleration  due  to  gravity  is  modeled  in  the  vector 
(gx,  9y»  9z^T*  This  gravity  model  includes  the  second 
harmonic  of  the  gravity  potential.  Since  the  coordinate  system 
is  moving  with  the  earth,  there  is  an  apparent  Coriolis 
acceleration  (Acx,  ACy,  ACZ)T  in  the  equation.  The  third 
term  in  each  acceleration  component  of  f  is  for  aerodynamic  drag 
on  the  target.  The  factor  D  in  the  drag  terms  is  used  to  model 
the  effects  of  atmospheric  pressure  and  the  ballistic  coefficient 


of  the  target.  It  is  computed  in  real  time  as  a  function  of 
altitude  from  a  table  of  quadratic  polynomials.  Each  polynomial 
is  used  in  a  specified  altitude  interval.  This  table  of 
polynomials  and  intervals  is  prepared  from  a  nominal  trajectory 
set  up  prior  to  a  mission. 

For  more  mathematical  details  on  these  equations  of  motion 
see  Appendix  A. 

2.3.  Measurements  Used  by  the  Filter 

The  measurements  used  by  the  filter  are  estimates  of  TRADEX 
range,  azimuth,  elevation,  the  Gellinam  bistatic  range  sum,  and 
the  Illeginni  bistatic  range  sum.  The  three  range  measurements 
are  assumed  to  have  been  calibrated  and  corrected  for  all 
significant  distinctions  between  the  apparent  propagation  times 
and  the  geometric  distances  between  the  target  and  the  three 
sites. 

Currently  the  real  time  program  corrects  the  range  m 
measurements  for  tropospheric  refraction  and  for  the  delay 
induced  by  range-doppler  coupling  in  the  pulse  compression 
network . 

The  time  of  validity  of  the  five  measurements  is  the 
estimated  time  of  incidence  of  the  radar  pulse  on  the  target. 
Since  TRADEX  measures  its  azimuth  and  elevation  at  a  slightly 
different  time,  the  real  time  program  extrapolates  these 


angle  measurements  to  the  time  of  incidents  using  the  angular 
rates  and  accelerations  reported  to  it  by  TRADEX. 

The  filter  adjusts  its  estimate  of  target  position  and 
velocity  on  the  basis  of  the  difference  between  the  calibrated 
and  corrected  measurements  described  above  and  corresponding 
values  computed  from  the  predicted  position.  This  predicted 
position  is  determined  from  previous  measurements.  The  formulas 
for  these  computed  values  are  given  in  Appendix  C. 

The  tristatic  filter  updates  its  estimate  of  the  target 
position  and  velocity  once  every  100  millisecond  cycle  of  the 
KDC  real  time  program.  During  each  cycle,  any  subset  of  the  five 
corrected  measurements  described  above  may  be  processed  by  the 
filter.  Typically  this  subset  consists  of  TRADEX  range  and  both 
bistatic  range  sums.  The  choice  of  subset  may  be  altered  for 
each  update  cycle.  This  capability  allows  the  real  time  program 
to  avoid  processing  bad  measurements  and  to  start  the  filter  with 
TRADEX  range  and  angle  measurements  before  the  remote  sites 
acquire  the  target. 


3.  ALGORITHMS  AND  EQUATIONS  USED  IN  THE  FILTER 

This  section  describes  the  algorithms  and  equations  used  to 
implement  the  tristatic  filter  based  on  the  mathematical  model 
described  in  Section  2.  The  next  subsection  gives  an  overview  of 
the  structure  of  the  filter. 

3.1.  Overview  of  the  Filter 

The  general  organization  of  the  tristatic  filter  is  shown  in 
Fig.  2.  During  each  100  millisecond  cycle  of  the  real  time 
program,  each  block  is  executed  once. 

The  Raw  Measurements  consist  of  the  TRADEX  range,  azimuth 
and  elevation  as  well  as  the  propagation  times  for  the  bistatic 
measurement  to  each  remote  site.  The  "Apply  Corrections"  block 
calibrates  the  Raw  Measurements  and  makes  the  adjustments 
discussed  in  Section  2.3. 

The  resulting  Corrected  Measurements  are  estimates  of  the 
geometric  TRADEX  range,  azimuth,  elevation,  Gellinam  bistatic 
range  sum,  and  Illeginr.i  bistatic  range  sum  at  time  of  target 
illumination.  These  five  corrected  estimates  are  termed 
geometric  ones  because  they  correspond  to  the  values  derived  from 
the  locations  of  the  sites  and  target  using  geometry. 

The  Prefiltet  block  in  the  diagram  checks  the  Corrected 
Measurements  for  reasonableness  by  comparing  them  with  the 
previous  measurements.  The  algorithm  for  this  comparison  is 
discussed  in  Section  3.2. 


The  Valid  Measurements  are  sent  on  to  the  Sequential  Update 
block  along  with  an  indication  of  which  measurements  have  been 
selected  for  use  in  the  update  process. 

Prior  to  executing  the  Sequential  Update  block,  the  filter 
predicts  the  target  position  and  velocity  (Predicted  State  in 
Fig.  2)  using  the  Previous  State  vector  derived  from  previous 
measurements.  This  is  done  by  integrating  the  equations  of 
motion  of  the  target  over  the  tiro,-  interval  between  the  time  of 
va  idity  of  the  previous  measurements  and  the  time  of  validity  of 
the  current  measurements.  The  accuracy  of  the  Predicted  State  is 
given  in  terms  of  the  Predicted  Covariance  Matrix,  which  is 
determined  from  the  Previous  Covariance  Matrix,  the  Process  Noise 
Covariance  0,  and  the  equations  of  motion.  The  initial 
conditions  for  this  integration  are  the  target  position  and 
velocity  that  constitute  the  Previous  State  in  Fig.  2.  The 
"Predict"  block  also  extrapolates  the  Previous  Covariance  Matrix 
to  determine  the  Covariance  Matrix  for  the  Predicted  State. 

The  Sequential  Update  combines  each  Valid  Measurement  with 
the  Predicted  State  vector  to  improve  the  accuracy  of  the 
target  position  and  velocity  estimates.  These  estimates 
constitute  the  Updated  State.  The  Updated  Covariance  Matrix 
indicates  the  accuracy  of  the  Updated  State.  Measurement 
Variances  are  used  to  control  the  weights  applied  to  the 
observation  to  obtain  the  Updated  State.  If  the  Prefilter 
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selects  no  measurements  for  use  in  the  Sequential  Update,  then 
the  Predicted  State  and  covariance  matrix  become  the  updated 
ones.  This  effects  a  coast  of  approximately  200  milliseconds  that 
might  be  extended  during  the  processing  of  the  next  measurement. 

To  get  started,  the  filter  needs  the  following  things 
supplied  to  it. 

(1)  Target  position  and  velocity  to  use  as  a  "Previous 
State"  vector. 

(2)  A  Covariance  Matrix  for  the  error  in  this  "Previous 
State"  vector. 

(3)  The  Process  Noise  Covariance  Matrix  Q. 

(4)  The  measurement  variances  used  in  the  Sequential 
Update, 

The  "Previous  State"  vector  is  taken  from  an  already  existing 
track  file,  such  as  the  TRADEX  monostatic  track  file.  The  other 
items  in  the  above  list  are  supplied  to  the  filter  by  operators 
when  the  real  time  operations  start, 

3.2.  The  Prefilter 

The  job  of  the  prefilter  is  to  detect  glitches  in  the  data 
so  that  only  valid  measurements  are  fed  to  the  Sequential 
Update.  The  difference  is  computed  between  each  actual  range 
measurement  and  the  corresponding  value  computed  from  the 
predicted  target  position  for  the  same  time.  The  average  and 
standard  deviation  are  computed  for  the  previous  30  such 


residuals,  and  the  residual  for  the  current  measurement  is 
compared  to  them.  A  detailed  statement  of  the  glitch-catching 
algorithm  for  TRADEX  range  follows. 

Let  tfc  be  the  time  of  validity  of  the  k-th  TRADEX  range  in 
the  sequence  of  range  values  being  processed. 

Let  Rp( k )  =  Rt  (?(tk  |  tfc-1)),  the  predicted  TRADEX 

range  computed  from  the  predicted  state  vector  If  (t^  |  tj^-i). 

Let  Rj^(k)  =  the  measured  TRADEX  range  at  time  t^. 

For  k  =  1  to  30  define 

d(k)  =  Rm  (k)  -  Rp  (k) 

For  each  k  >  31  define  the  average  d(k)  and  standard 
deviation  a(k)  of  the  previous  30  residuals  by 


If 


d(k)  =  1 


30 


c2(k)  =  ^ 


RM  (k)  - 


30 

Z 

1=1 

30 

Z 

1=1 


Rp  (k) 


d(k  -  *) 

(d(k  -  i)  -  d  (k))2 
-  d  (k)  |  <  5  a ( k ) 


then 

(1)  Use  R^  (k)  in  Sequential  Update  for  time  t^ 
(i.e.  Rm  (k)  is  Valid  Measurement) 

( 2 )  Define  d ( k )  =  R^  { k )  -  Rp  ( k ) 


Otherwise 
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(1) 


Declare  a  glitch  and  do  not  use  (k)  in  Sequential 
Update 
(2)  Define  d(k)  by 


( 

-500 

i 

if 

rM 

( k )  —  Rp 

(k) 

< 

-500 

d(k)  =  ! 

\  RM  (k)  -  Rp  (k) 

if 

rM 

(k)  -  Rp 

(k) 

<_ 

500 

l  500 

if 

rM 

(k)  -  Rp 

(k) 

> 

500 

The  above  algorithm  is  written  for  TRADEX  range. 

To  monitor  the  Gellinam  bistatic  range  sum  measurement,  a 
similiar  algorithm  is  used.  In  this  version  the  measured 
bistatic  range  sum  replaces  the  measured  TRADEX  range  for  each 
time  tk,  and  Sq  (1c  (tk  |  i ) )  replaces 

Rt  (xV-k  j  tk-i)).  A  similiar  algorithm  is  also  used  for 
the  Illeginni  bistatic  range  sum. 

3.3.  Predicting  Target  Position  and  Velocity 
The  prediction  of  the  target  position  and  its  velocity  are 
based  on  numerical  integration  of  the  equations  of  motion 
described  in  Section  2.  The  numerical  integration  consists  of 
one  step  of  the  Runge-Kutta  method  from  the  time  of  the  previous 
measurement  tk_i  to  the  time  of  the  current  measurement  tk. 

The  initial  condition  for  the  integration  is  the  result  of  the 
previous  Sequential  Update,  namely  x  (tk-i  tk-i). 

The  estimate  of  the  state  vector  at  time  tk  is 
x  (fck  |  fck- 1 )  =  x  (tk-i  |  tk_i)  +  (f\  +  2 Tz  +  2?3  +t)/6 
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where 

?!  =  At  f  (  X  ( tk_ !  |  tk_i)) 

C2  s  At  f  (  X  (tk_x  tk_i)  +  1/2  ?i) 

?3  =  At  f  (  x  (tk_i  tk_x)  +  1/2  ?2 ) 

>  ^  ^  .  ■$> 

Ci*  »  At  f  (  x  (tk_i  I  tk_ x )  +  ? 3  ) 

3.4.  Prediction  of  the  Covariance  Matrix 
The  computation  of  the  covariance  matrix  associated  with 

I  i 

x  (tk  tk_i),  namely  P  (tk  tk_x)  i-s  based  on  the 
equation 

P  <fck  |  fck-l)  *  *(tk)  P  <  tk_ i  |  tk_  x )  *(tk)T 
+  *(tk)  0  *(tk)T. 

In  this  equation  #(tk)  is  an  approximate  transition  matrix 
given  by 

•$> 

*(tk)  =  I  +  AtF  (X(tk  I  tk_x>) 

> 

^  3f  > 

where  F  (x)  =  -3-  (x), 

oX 

is  the  Jacobian  of  the  vector  function  in  the  right  hand  side  of 
the  equations  of  motion.  Using  this  formula  for  the  transition 
matrix,  we  can  express  the  predicted  covariance  matrix  as 
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p  <tk 
where 


|  fck-l)  *  P  +  At(F  P  +  (FP)t  +  At  FPFt) 

P  =  P<tk  |  fck>  +  Of 

F  =  F(x  (tk  |  tk)), 
and  At  =  tk  -  tk_j. 

The  filter  currently  uses  this  formula  for 
P  (tk  j  tk_i)  with  a  simplified  version  of  the  Jacobian  F, 
instead  of  the  exact  Jacobian  matrix. 

3.5.  Using  the  Simplified  Jacobian  to  Predict'  the 
Covariance  Matrix 

This  section  discusses  the  simplifications  applied  to  the 

df 

Jacobian  matrix  -3.  and  the  algorithm  that  takes  advantage  of  the 

3x 

them. 

One  simplification  is  to  replace  the  nine  partial 

derivatives  of  f4 ,  f5 ,  and  f6  with  respect  to  x,  y,  and  z  in  the 

Jacobian  matrix  with  zeroes.  Another  simplication  occurs  in  the 

expressions  for  — — — -y  ! — -5--.  and  ■■■  —  &- 

3y  3z  3x  3x 

where  the  rotation  rate  of  the  earth  u>  is  taken  to  be  zero. 

These  simplications  eliminate  all  the  terms  in  the  partial 
derivatives  of  acceleration  except  those  involving  drag. 

The  simplified  Jacobian  can  be  written 
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where  each  block  matrix  is  3  by  3  and  the  drag  matrix  A  is  given 
by 


The  algorithm  reduces  the  amount  of  computation  required  to 
predict  the  covariance  matrix  by  taking  advantage  of  the  blocks 
of  zeros  and  the  identity  matrix  in  the  simplified  Jacobian 
matrix  as  described  below. 

Partition  the  covariance  matrix  P  into  3  by  3  blocks. 


where  B  and  E  are  symmetric  matrices. 

The  needed  matrix  products  are  given  by: 


Let  M  =  ACT,  n  =  AE  and  K  =  AEAT.  Then  we  have 


The  algorithm  consists  of  the  following  steps: 

(1)  Compute  the  drag  term  A  using  equation  2. 

(2)  Compute  the  3  by  6  matrix  [  M  N  ]  using  equation  6. 

(3)  Compute  K  =  N  AT. 


(4)  Construct  the  matrix  products  using  equations  7,  8  and  9 
and  the  entries  in  M,  N  and  K. 


(5)  Compute  the  Predicted  Covariance  Matrix  P(tjc  |  t^..^), 

from  last  equation  in  previous  section. 

3.6.  The  Sequential  Update 

The  Sequential  Update  combines  the  MMS  measurements  and  the 
predicted  target  position  and  velocity  to  improve  the  estimates 
of  position  and  velocity. 

The  Sequential  Update  consists  of  a  sequence  of  Kalman 
updates,  each  of  which  combines  a  single  scalar  measurement  with 
target  position  and  velocity  estimates.  This  approach  presumes 
that  the  MMS  measurements  at  any  one  time  are  uncorrelated. 

The  measurement  selection  in  the  prefilter  determines  which 
scalar  updates  are  done  during  a  100  millisecond  cycle  of  the 
real  time  program.  Fig.  3  shows  all  five  scalar  updates  that  are 
possible  and  their  inter-connections.  Fig.  3  assumes  that  all 
five  measurements  have  been  selected  for  processing.  More 
typically  only  the  three  MMS  measurements  are  selected,  in  which 
case  the  scalar  updates  for  TRADEX  azimuth  and  elevation  are 
deleted  from  the  Sequential  Update  as  shown  in  Fig.  4.  Another 
useful  configuration  for  the  sequential  update  is  one  which 
processes  only  TRADEX  data.  This  is  shown  in  Fig.  5.  Such  a 
configuration  could  be  used  to  start  the  filter  before  the  remote 
sites  acquire  the  target. 

The  configurations  for  the  Sequential  Update  shown  in  the 
figures  presume  that  the  pretilter  has  not  detected  a  glitch  in 
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Fig.  3.  Sequential  update  using  all  five  measurements. 
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SG 
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Fig.  4.  Sequential  update  using  MMS  data. 


Fig.  5.  Sequential  update  using  TRADEX  data 


Fig.  6.  Generic  form  of  scalar  update. 


the  MMS  data.  If  a  glitch  is  detected  by  the  Prefilter,  the 
measurement  selection  passed  to  the  sequential  update  block  in 
Fig.  2  will  exclude  the  bad  scalar  measurement.  In  reaction  tc 
this  deletion,  the  sequential  update  will  also  exclude  the 
corresponding  scalar  update  from  the  configuration  chosen  a 
priori.  Unless  the  Prefilter  detects  a  second  glitch  in  the 
following  100  millisecond  cycle,  the  associated  scalar  update  is 
put  back  into  the  sequence  of  scalar  updates. 

Each  box  in  Figs.  3,  4,  and  5  represents  a  scalar  update 
based  on  the  Kalman  update  equations.  A  generic  version  of  this 
scalar  update  is  shown  in  Fig.  6.  In  terms  of  the  notation  in 
Fig.  6,  the  scalar  update  equations  are: 

y j+i  =  yj  +  Kj  -  fcj  (yj)) 

Uj+1  =  (I  -  K  Hj)  Uj 
Kj  =  Uj  Hj  (Hj  Uj  HjT  +  Oj2)-1 
Here  yj  is  the  state  vector  detemined  from  previous 
measurements,  nij  is  the  scalar  measurement,  hj  is  the 
associated  measurement  functicn,  Uj  the  covariance  matrix  for 

yj,  Hj  is  the  Jacobian  of  h j ,  namely 

ahi  > 

H-;  =  i  evaluated  at  x  =  Vi 

J  3x  J 

and  oj  is  the  standard  deviation  associated  with  the 

The  matrix  Kj  is  the  gain  matrix.  In  Fig.  3 


measurement  mj 


^  <fck  |  tk_!)  is  to,  P  <tk  |  tk_i)  is  uo 

t  (tk  |  tk)  is  y5  and  P  (tk  |  tk)  isus* 

Similiar  associations  hold  for  other  configurations  of  the 
Sequential  Update. 
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4. 


EXAMPLE  OF  TRACKING  WIT*  TRISTATIC  FILTER 


This  example  is  drawn  from  the  mission  GT100.  The  iriscatic 
filter  estimates  of  target  position  and  velocity  can  be  used  to 
calculate  the  strobe  times  for  the  remote  sites.  Estimates 
generated  by  the  TRADEX  monostatic  tracker  can  also  be  used  for 
this  purpose.  The  difference  between  the  expected  time  of 
arrival  of  the  peak  return,  as  determined  by  the  tracking  filter, 
and  the  actual  time  for  it  are  plotted  for  each  site  in  Figs.  7, 
8,  and  9.  Prior  to  1662  seconds  after  launch  the  TRADEX 
monostatic  tracker  is  directing  the  laying  of  the  strobes  at  the 
remote  sites.  After  the  1662  second  mark,  the  tristatic  filter 
is  directing  the  remote  sites.  The  strobes  at  TRADEX  are  always 
directed  by  the  TRADEX  monostatic  tracker.  These  plots  show  a 
dramatic  improvement  in  the  accuracy  of  the  laying  of  strobes 
when  controlled  by  the  tristatic  filter. 


TIME  AFTER  LAUNCH  (s) 


Fig.  7.  Strobe  laying  error  for  TRADEX  range. 
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APPENDIX  A  DETAILS  OF  THE  EQUATIONS  OF  MOTION 

In  Section  2  the  equations  of  motion  of  the  target  are 
discussed.  This  appendix  provides  the  mathematical  details.  The 
acceleration  due  to  gravity  acting  on  the  target  at  {/.,  y,  z)  is 
(Qx'  Qy>  9z )T  where 
9x  =  x  9c 

9y  =  [y  -  Rs  sin  <3  -  a)J  gc  +  cos  3  g^ 

g?  =  [  z  +  H  +  Rs  cos  (3  -  a ) ]  gc  +  sin  3  g<j 

9c  =  “T  [1  -  J2  (  ^  )2  P3'(sin  Y)] 

9d  =  J2  (  ^  )2  p2  ’  <sin  Y  ) 

R  =  distance  of  target  from  center  of  earth 
R2  =  x2  +  [y  -  Rs  sin  (3  -  a)]2  + 

[  z  +  H  +  Rg  cos  (3  “  a) ]2 
Rg  =  Equatorial  radius  of  earth 
Rg  =  Radius  of  earth  at  TRADEX 
H  =  height  of  TRADEX  above  earth 
u  =  Gravitational  constant  for  earth 
(3.986  x  lO20  cm3/sec2) 

Y  =  Geocentric  latitude  of  target 

sin  y  =  [ Rs  sin  y  +  H  sin  3  +  y  cos  3  +  z  sin  3]  /  R 
J2  =  Second  harmonic  of  earth's  gravity  potential 

Y  =  Geocentric  latitude  of  TRADEX 

3  =  Geodetic  latitude  of  TRADEX 
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P2  and  P3  are  the  second  and  third  order  Legendre 
polynomials.  Their  derivatives  are: 

P2'(X)=  3X,  P3'(X)  =  1/2  (15  X2  -  3) 

The  target's  motion  and  the  earth's  rotation  lead  to  the 
Coriolis  acceleration  components, 

Acx  =  -  2  a)  z  cos  M  2  u  y  sin  0  + 

ACy  =  -  2  a)  x  sin  3  +  ui2y  sin2  3 

-  oj2  sin  $  [R  cos  a  +  (H  +  z)  cos  3] 

AC2  =  2  (ox  cos  3  -  a)2y  sin  3  cos  3 

+  u)2  cos  3  [R  cos  a  +  (H  +  z)  cos  3] 

where 

ai  =  angular  velocity  of  the  earth. 

The  drag  model  is  explained  in  Section  2 . 
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APPENDIX  B 


JACOBIAN  OF  EQUATIONS  OF  MOTION 


The  Jacobian  of  the  right  hand  side  of  the  system  of 


Hr 

equations  of  motion  is  ( 
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APPENDIX  C 


MEASUREMENTS  AS  FUNCTIONS  OF  TARGET  POSITION 


The  formulae  for  the  measurements  in  the  notation  in  Fig. 
1  are  the  following 
TRADEX  range  is 

RT  =  /  x2  +  y2  +  z2' 

TRADEX  azimuth  is 

A  =  arctan  (x/y) 

TRADEX  elevation  is 

z 

F  =  arctan 

/  x2  +  y2' 

The  Gellinam  bistatic  range  sum  is 


Sq  =  Rt  +  ^G 


where 

Rg  =  /  (x-xG)2  +  (y-yc)2  +  (z-zG)2k 
The  llleginni  bistatic  range  sum  is 


Si  =  R^  +  Rj 

where 

Rj  =  /  (x-xj;)2  +  (y-yj)2  +  (z-zj)2  ' 

The  Sequential  Update  in  the  tristatic  filter  uses  the 
partial  derivatives  of  these  five  measurement  functions.  These 
derivatives  are: 
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